Hybrid methods for witnessing entanglement in a microscopic-macroscopic system 
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We propose a hybrid approach to the experimental assessment of the genuine quantum features of a general 
system consisting of microscopic and macroscopic parts. We infer entanglement by combining dichotomic 
measurements on a bidimensional system and phase-space inference through the Wigner distribution associated 
with the macroscopic component of the state. As a benchmark, we investigate the feasibility of our proposal in 
a bipartite-entangled state composed of a single-photon and a multiphoton field. Our analysis shows that, under 
ideal conditions, maximal violation of a Clauser-Horne-Shimony-Holt-based inequality is achievable regardless 
of the number of photons in the macroscopic part of the state. The difficulty in observing entanglement when 
losses and detection inefficiency are included can be overcome by using a hybrid entanglement witness that 
allows efficient correction for losses in the few-photon regime. 



I. INTRODUCTION 

An open challenge for fundamental quantum physics is to 
affirm the quantum nature of a system that puts together a 
microscopic part and a mesoscopic one. This hybrid sce- 
nario can emerge in completely different experimental plat- 
forms ranging from individual spin systems interacting with 
multimode cavity fields, such as transmon qubits in coplanar 
transmission-line resonators 111 01 ^ to ionic impurities embed- 
ded in ultracold atomic samples, such as the systems con- 
sidered in some recent experiments reported in iH 0] ■ An- 
other possible physical approach exploits a massive tiny mir- 
ror interacting optomechanically with a single photon within 
a Michelson interferometer ||5|-l21 • This endeavor could con- 
tribute to challenge the observability of quantum features at 
the macroscopic level, which is one of the most fascinat- 
ing open problems in quantum physics. The difficulties in- 
herent in such a quest are manifold, and they are related 
on the one hand to the unavoidable interaction of the sys- 
tem with the surrounding environment ifiol - fisll . On the other 
hand, one faces the debated problem of achieving a measure- 
ment precision sufficient to observe quantum effects at such 
macroscales ifTilfTsll . In this context, it has been experimen- 
tally proven that a dichotomic measurement performed upon 
a multiphoton-entangled state is not sufficient to catch quan- 
tumness lfl6ll . The accuracy of the measurement is crucial for 
the observation of quantum features and should be put on the 
same footing as the use of proper entanglement and nonlocal- 
ity criteria for macroscopic quantum systems lfT5l[T7l - l2lll . 

To successfully tackle the manipulation and characteriza- 
tion of hybrid systems the following question is still open: 
How can we ascertain the nonclassical nature of a multipar- 
tite state that, per se, does not meet the criteria for quantum- 
ness that have been designed for system components of equal 
dimensionality? Our work provides a quantitative answer to 
this broad question. We introduce an investigative platform 
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that can be built up without the necessity for information on 
the state itself, and this supports the general validity and broad 
applicability of our results. 

We introduce a hybrid method to demonstrate experimen- 
tally the truly quantum mechanical features of a general 
microscopic-macroscopic system beyond any assumption on 
its state and without the necessity of any a priori state knowl- 
edge. We infer the entanglement properties by means of a 
hybrid approach that combines dichotomic measurements on 
a bidimensional system and phase-space inferences through 
the Wigner distribution associated with the macroscopic com- 
ponent of the state. Here, through the use of a hybrid en- 
tanglement test, we identify a valuable tool for our goals. 
While the microscopic part of the state is measured using spin- 
1/2 projection operators, the macroscopic counterpart under- 
goes phase-space measurements based on the properties of 
its Wigner function iITtIi . At variance with previous propos- 
als lfl7il22ll . the approach presented in this paper is tailored to 
fully exploit the polarization-spin degree of freedom on both 
the microscopic and the macroscopic subsystems. We ana- 
lyze the effects of losses on a Clauser-Horne-Shimony-Holt- 
like (CHSH-Uke) inequality test ll23ll and show that maximum 
violation is achieved when losses are absent, regardless of the 
size of the macroscopic part of the state. This is not the case 
under nonideal conditions. However, we show how losses can 
be efficiently taken into account so as to infer entanglement of 
our multiphoton state. 

As a paradigmatic microscopic-macroscopic system 
(MMS), we investigate the state obtained from a fully 
microscopic-entangled system through an amplification 
process 11241 Izsll . Such system has been further consid- 
ered recently as a benchmark to perform nonlocality tests 
with human-eye threshold detectors ll26ll or as a platform 
for absolute radiometry IztIi . At variance with respect to 
Refs. ll24l[25ll . our approach does not require any assumption 
on the system under investigation and hence represents a 
genuine entanglement test. 

The present paper is organized as follows. In Sec. HIl we 
introduce and define the CHSH-based entanglement inequal- 
ity based on hybrid measurements for the single-photon and 
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the multiphoton modes. Then, in Sec. |III]we discuss how the 
CHSH-based test defined in Sec. |II]can be modified to ob- 
tain an entanglement witness tailored to be applied in a lossy 
scenario. Finally, in Sec. |IV]we provide a specific example 
of a joint optical system composed by a single-photon and a 
multiphoton field based on the process of optical parametric 
amplification. We then run both the CHSH-based test and the 
entanglement witness on this system in order to identify in 
which range of the system's parameters the entanglement can 
be addressed with our approach. 



II. HYBRID ENTANGLEMENT TEST BASED ON BELL'S 
INEQUALITIES 

Let us consider a general MMS state with its microscopic 
part embodied by a single-photon polarization state (a qubit). 
We take the macroscopic part, on the other hand, as encoded 
in the multiphoton state of a continuous-variable (CV) sys- 
tem. The two subsystems are supposed to be entangled by a 
mechanism whose details are inessential for our tasks here. 
A benchmark state of such situation will be provided later 
Polarization measurements performed over the state of the 
single-photon mode are described by the Pauli spin op- 
erator a-^icj)) = \(f)) a{4'\ ~ \4'i)a{4'i^V where is the direc- 
tion identifying the polarization state in the Poincare sphere 
and is its orthogonal direction. The CV measurements, on 
the other hand, are given by i^'^^^ (ctx' x)~^x i'^x)®^x±' 
where tif {a,)^bf {a^){~lf^ Df\a^) is the displaced 
parity operator built from the displacement Df{ai) (ctiGC) 
and the number operator hf (i={Xi X-l} stands for the polar- 
ization state. Such operators can be directly measured ll28ll29ll 
by combining the input field with a coherent state in a low re- 
flectivity beam splitter and by measuring the parity of the out- 
put field [see Fig. [11(a)]. However, such technique requires a 
photon-counting technique with very high efficiency, a condi- 
tion extremely difficult to achieve with the present technology. 

An indirect measurement of the average value of the 
displaced parity operators can be performed by exploiting 
the connection between (n(a)) and the Wigner function of 
the state [see Fig. [T] (b)]. Indeed, the average value of 
the measurement operator on state pf of the multiphoton 
mode is related to the value of its Wigner function at a^, 

M^|(a,)=(2/7r)Tr[nf (aOpf ]■ The latter can be easily re- 
constructed using homodyne measurements. We define the 
qubit-CV coiTclator C(ax, %; (/))=((T'^(0)®n^,^^ (a^, x)), 
which is evaluated on a general MMS pab, and the CHSH- 
based entanglement parameter 

B=c{a'^,^- x'; 0)+CK, x; 0')-CK. x; <t>)- 

(1) 

A more detail discussion can be found in Appendix [A] 
As the average values of the outcomes of the (T"^(0) and 
xxi('^X'X) measurements is limited by {<7^{(f))) < 1 e 
(fi^^^ (a^, x)> < 1, for ah separable states flie bound \B\<2 
holds. A violation of this bound witnesses an entangled state. 
The measurement settings for the single-photon mode k^i 
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FIG. 1. (Color online) Hybrid entanglement test on an optical MMS 
state generated by a "black-box." The single-photon mode is 
measured by a polarization-detection apparatus, while the multipho- 
ton mode ks undergoes polarization projection and the measure- 
ment of the displaced parity operators, (a) Direct measurement of 
the 11(a) displaced parity operators, (b) Indirect measurement of the 
average value (n(a)) of the displaced parity operators by exploiting 
a homodyne-detection apparatus. 



[multiphoton mode ks] are given by the measured polariza- 
tions ((/), (j)') [measured polarizations (x, x') and the chosen 
phase-space points (a^, ap]. This requires a standard polar- 
ization detection system for the microscopic mode and a ho- 
modyne detection system for the multiphoton one, as shown 
in the scheme presented in Fig.[T](b). 

We conclude by observing that the inequality of Eq.([T]i be- 
comes a nonlocality test when the displaced parity operator 
are directly measured on the multiphoton field [see Fig. [T](a)], 
since no assumption are necessary on the detection appara- 
tus. In this case, the outcome of the o'"^(0) and H^-^^ i'^x^ x) 
measurements can only be ±1, and the use of a local-hidden- 
variable (LHV) model imposes the bound |/B|<2 |[23il on the 
B parameter A violation of this bound confutes all LHV the- 
ories. 



III. HYBRID ENTANGLEMENT WITNESS WITH LOSSES 

The test presented above can be modified so as to embody 
a witness able to reveal entanglement when the state at hand 
is affected by losses. This is modeled by inserting a beam 
splitter of transmittivity ?7G[0, 1] in the path of the modes at 
hand, "tapping" the corresponding signal ifisll . The choice 
77=1 (7y=0) corresponds to a lossless (fully lossy) process. To 
this end, the measurement performed on the tt^ polarization 
of the multiphoton part is replaced by the operator 1I21I1 

o.^K;.)-0"^'^"^^+('"')'^ (2) 

i2nf(ax)-if if 776(0, 0.5]. 
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In this way, the overall measurement on the macroscopic 

subsystem reads d^^^^{a^,XiV)=0^io:x'^'n)'^^x±- 
For any separable state being measured after the lossy 

process, 1(0^,^^ X; ??)>„l < 1 S llH IM EH]. 
Hence, by introducing the MMS correlator 

Cr,(ax' X, 0)=('5-^(^)«'O^,x_L ("x' X; v))v^ we define 

Any separable state undergoing a lossy process on mode 
is bound to satisfy |yV^^''|<2 (see Appendix iBb. Violation of 
this inequality witnesses entanglement in the system. Such a 
bound can be explained by considering that separable states 
do not violate CHSH inequalities, and local processes such as 
losses cannot increase their nonlocal character. It is important 
to notice that, by virtue of the assumption that the macrostate 
of mode ks undergoes losses ij before (rather than at) de- 
tection, this entanglement witness reveals the presence of en- 
tanglement without any assumption on the MMS source (see 
FigHJ. On the other hand, the lossy mechanism can be shifted 
to occur just before measurement, thus modeling the effects of 
a non-ideal detector. For rj=l, W,, coincides with the CHSH- 
based parameter B in Eq. ([T]i. 



IV. EXPERIMENTAL BENCHMARK 

In this section, we analyze in details a specific optical sys- 
tem to evaluate the effectiveness of our hybrid approach. As 
a benchmark for the hybrid CHSH-based entanglement test 
and entanglement witness described above, we analyze the 
MMS-state-source addressed in Ref. ll24ll . A layout of the 
system is reported in Fig. |2] The polarization singlet state 
\^-)abH\H)a\V)b-\V)a\H)b)/V2 of a photon pair is 
generated in a nonlinear crystal through a spontaneous para- 
metric down-conversion (SPDC) process. Here \H) i\V)) 
stands for the horizontal (vertical) polarization state. The 
photon populating mode kg is then injected into an opti- 
cal parametric amplifier (OPA) in a collinear configuration. 
Since the OPA implements a unitary operation, the symme- 
try of \iP~)ab is preserved by the amplification process and 
the overall state |4'-)AB-(|0>^|$t>B-|0±)A|$'^)B)/V2 



polarization analysis 




FIG. 2. (Color online) Layout of the MMS source based on the pro- 
cess of optical parametric amplification of a single photon belonging 
to an entangled pair. 



maintains rotational invariance form for any polarization ba- 
sis. Here, 1$"^) are the multiphoton states generated by am- 
plification of a single -photon polarization state Quan- 
tum entanglement between the micropart and the macropart 
of \'i'^)AB has been demonstrated |l24il under a supplemen- 
tary assumption on the source 12511 . The OPA performs 
the optimal cloning process only for equatorial polarization 
7f0=(7r/f+e"^7r\/)/\/2- We thus restrict our attention to this 
subset of polarization states, which motivates our choice for 
(7"^(0) performed above. 




{1-V){n) 

FIG. 3. (Color online) CHSH-based parameter Bn as a function of 
the number of lost photons ( 1 — 77) (n) for different values of the gain 
g. We show the local realistic boundary Slhv=2. 

We now discuss the results of the CHSH-based test and 
the application of the entanglement witness to the MMS state 
given in Fig.|2] We begin analyzing the CHSH-based inequal- 
ity ([1]) in the lossless case (rj—l). The correlation operator 
evaluated on \'^~)ab takes the form (see Appendix |C]i 

ax^, P^,x; 0)-(l-2) cos[2(x-0)]e-^, (4) 

where Z=2(e~^^X^+e^ffP^) is a function of the ro- 
tated variables X^=X;^ cos(x/2)— sin(x/2) and 
Px=^x sin(x/2)+Px cos(x/2). {X^,P^) are the field 
quadratures and ay^=X^+iP^. The correlator in Eq. (|4| 
is maximized at the origin of the phase space, where 
C(0, 0, x; 9^)= cos[2(x — </>)], which is independent of the 
gain of the amplifier g and the number of generated pho- 
tons n — sinli'^ g. The correlator has the same form as a 




g g 



FIG. 4. (Color online) (a) Contour plot of the shifted loss function 
£{g,r])—2^^^^ as a function of the gain g and the detection effi- 
ciency 77. (b) Contour plot of the negativity of the Wigner function 
of an amplified single-photon state 1331 against g and 77, evaluated at 
the origin of the phase space. In both panels the solid line divides 
the region of entanglement (|S,,|>2, above the line) from the one in 
which entanglement cannot be inferred (|S,, | < 2, below the line). 
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CHSH-Based Entanglement test: YES 
Entanglement Witness: YES 
CHSH-Based Entanglement test: NO] 
Entanglement Witness: YES 
CHSH-Based Entanglement test: NOl 
Entanglcntcnt Witness: NO 



FIG. 5. (Color online) (a) Plot of Wn vs the detection efficiency rj and the nonlinear gain g. (b) Contour plot of the effective loss function 
h(ri)C{'q, g). Entanglement can be revealed in the region above the black line, (c) Summary of the results obtained from our tests. We identify 
three regions in the [rj, g) space, depending on whether entanglement can be demonstrated with our techniques. 



Bell-CHSH test performed on a polarization photon pair, 
where spin- 1/2 operators are measured. The CHSH-based 
parameter B is then maximized by choosing the measurement 
settings for (0, 0, x, x') corresponding to such case, which 
ensures the maximum degree of violation of the local realistic 
boundary, i.e., B=2\/2. 

We are now in a position to address the possibility to ob- 
serve MMS entanglement under realistic experimental condi- 
tions. We thus analyze the effects of detection efficiency at 
the homodyne apparatus, while other sources of experimen- 
tal imperfections, as well as a more detailed derivation, are 
discussed in Appendix IdIi. The measurement of the gener- 
alized parity operator on the multiphoton mode can be 
performed using homodyne detection. Furthermore we in- 
clude in the qubit-CV correlator the possibility of a nonuni- 
tary detection efficiency on mode (see Appendix|D]i. By 
restricting our attention to the origin of the phase space, 
where maximum nonclassical effects are achieved, we get 
C,(0, 0, x; 4>)= cos[2(x-0)]/:(.g, v\ where 



7/[l + 2n(l - 7?)] 
(1 + 47?(1 - ?7)rT)3/2 



(5) 



is a loss function for the test. Hence, the maximum amount 
of violation is directly determined by the loss function as 
Bn=BC[g, rj). In Fig. [3] we show the value of Brj as a func- 
tion of the average number of lost photons, {l—rf){n), where 
(n) = 3?T+ 1 is the mean number of the generated photons af- 
ter the amplification process. The CHSH-based inequality of 
Eq. ([TJ is satisfied when only a moderate number of photons is 
lost. A lower bound 77/i,„=l/V2 for the detection efficiency 
can be found below, where a violation is no longer observed. 
On the other hand, at set values of i] there is a minimum gain 
giimiv) above which the presented test cannot detect MMS- 
entangled correlations. Such threshold value decreases with 
the reduction of the efficiency ?y. The behavior of in the 
{i], g)-plane is shown by the contour plot in Fig.|4](a). In or- 
der to relate the violation of the CHSH-based inequality to 
intrinsically nonclassical features enforced at the level of the 
macro-part of the state. Fig. |4](b) reports the negativity of the 
Wigner function of an amplified single-photon state versus rj 
and g [32]. We observe that the transition of B,^ to the region 
below the classical limit is directly linked to the decrease in 
the negativity of the Wigner function itself. Indeed the value 



of the MMS correlator is determined by the excursion of 
the Wigner function in X^=P^=0, as a function of the polar- 
ization of the injected photon. 

We complement the analysis of our MMS by discussing the 
use of the entanglement witness described above. The evalua- 
tion of the correlation operator over state |5'^)yi_B after losses 
leads to C,^{a^,x,(t))^Hv)C,,ia^,x;(t)), where h{ri)=l/-q 
(h{r])=2) for 1/2<77<1 (0<77<l/2). More details can be 
found in Appendix |E] Therefore, the entanglement witness 
can be directly obtained from the CHSH-based parameter as 
Wri~h{rj)Bri- InFig.|5](a) we report the dependence of as 
a function of 7] and g. For single-photon states (i.e., atg = 0), 
the correction of losses introduced by the factor h{r]) allows 
one to observe MMS entanglement up to As the num- 

ber of photons in the macrostate increases, the damping in 
the negativity of the Wigner function induced by losses scales 
more rapidly than 77 and the /i (77) -correcting term becomes 
less effective. Fig.|5](b) shows the behavior of the effective 
overall loss function h{ri)C{ri, g), highlighting the thresholds 
in g and 77, above which entanglement is observed. We note 
the non-monotonic behaviour obtained for the inefficiency pa- 
rameter at 77 = 0.5, which is a property of the witness itself. 
However, being Eq. (|2]i a witness for entanglement, no spe- 
cial meaning can be attached to the lack of violation of the 
separability condition jW^'^^l < 2. 



V. CONCLUSIONS AND PERSPECTIVES 

We have proposed an experimentally oriented approach to 
detect entanglement in a MMS -entangled state involving a 
single-photon and a multiphoton bipartite system. We have 
used a hybrid CHSH-based inequality and an entanglement 
witness, whose use against such a class of states is effective. 
Furthermore, the CHSH-based inequality can be adopted as a 
genuine nonlocality test when a direct measurement of the dis- 
placed parity operators is performed on the multiphoton field. 
As an experimental benchmark, we applied the proposed in- 
equalities to the bipartite state obtained by amplification of an 
entangled single-photon-singlet state. While our study spurs 
further interest in the identification of suitable tests in the 
high-loss and large-photon-number region, it paves the way 
to an experimentally feasible demonstration of entanglement 
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properties in an interesting class of states lying at the very 
border between quantum and classical domains. 



Appendix B: Hybrid polarization-continuous variables 
entanglement witness with inefficient detectors 
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Appendix A: Hybrid polarization-continuous variables 
CHSH-based test 



In this section we review the CHSH-based inequality per- 
formed in the paper. Our test is the extension of the Bell's 
inequality test proposed by Wodkiewicz in Ref. [17]. We be- 
gin by focusing our attention on the multiphoton mode k^. 
Our MMS, which is generated by amplification of an entan- 
gled polarization photon pair, is strongly correlated in such a 
degree of freedom. To exploit it, we define the measurement 
operator of the multiphoton state as 



(Al) 



Here nf = Df {a^){-lf^ Df\a^) is the generalized 
parity operator, where D^{ai) is the displacement operator 
and the subscript i = {x,X-l} describes the polarization 
mode. This definition of the measurement operator corre- 
sponds to the application of a displacement operator Df{ai) 
followed by a parity measurement. 

In order to detect the correlations present in the system in 
the polarization degree of freedom, we perform a measure- 
ment of the Pauli operator a^{4') on a single-photon mode. 
Here, o"((/)) is the Pauli operator along the direction of 
the Bloch sphere identified by the equatorial polarization state 
7?^ = 2~^/^(7f// + e"*7fy). The correlation of the joint system 
is then defined as 



(A2) 



where the averages are evaluated on the investigated \^)ab 
MMS state. Since this correlation operator corresponds to a 
set of dichotomic measurements, we can use the CHSH-based 
inequality [23] 

B=c(a;^, x'; </'')+cK, x'; 0)+c(ax, x; x; </'), 

(A3) 

Here, the measurement settings for the single-photon mode 
are given by the measured polarizations (0, 0'), while the 
measurement settings for the multiphoton mode k^ are given 
by the measured polarizations (x, x ) ™d the chosen phase 
space points (a^, a'^). 



In this section we discuss in details the hybrid entangle- 
ment witness defined in the paper Such an inequality is an 
extension of the CHSH-based test of Eq. (IA3b where differ- 
ent measurement operators are exploited in the multiphoton 
mode. The main idea of this extension is to take into account 
detection losses in order to build measurement operators apt 
for witnessing entanglement with an inefficient detection ap- 
paratus. To this end, the measurement performed on the 7?^ 
polarization of the multiphoton field can be replaced by the 
operator ll20l|2lll 



in^K) + (i-i)i^ ifi<^<i, 

2n^(aJ-i^ ifr;<i, 



(Bl) 

where rj is the detection efficiency of the apparatus. Such def- 
inition of the measurement operator is performed in order to 
correct the detrimental effect of losses on the properties of the 
detected state. Let us consider a general state |^)^ on spa- 
tial mode ks and polarization tt^. (Although we illustrate our 
argument using pure states of mode B, our arguments apply 
equally to mixed states). After losses occur, the state evolves 
into a density matrix pl^^. The average value of 0^{a^]iii) 
on such a density matrix gives 1I2TI1 



q-Wl^{a^)+[l-l] if i<r?<l, 



if < ry < i. 

(B2) 

Here, Wl^[pL^)\i, the Wi gner function of the detected state, 
which is related to the Wigner function of the initial state be- 
fore losses 1^")^ by the Gaussian convolution 



7r(l - 77) 



CX3 COO 



00 •/ —00 



X WiiX'P')e 



(B3) 

The measured Wigner function given in Eq. ( IB3b corre- 
sponds to the s-parametrized quasi-probability distribution 
Wi,{a^,s), of 1$)^ with s = Exploit- 
ing the properties of such distributions, it is straightforward to 
prove that 1I21I1 



(B4) 



for all values of rj. We can then define the overall measure- 
ment performed on the multiphoton state as 



Ox.xxK:x;'7) = o^K;'7)® i> 



(B5) 



with average values bounded by |(0^^^ (ttp,-, x; ?7))r;| < 1- 
The two-mode correlation operator for the entanglement wit- 
ness is then defined as: 



C'ia^.x; 0; v) = CT^(0) ® o^,x^ ("x. x; v) 



(B6) 
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where ct^ {(f)) is the PauH operator for mode k/i along the di- 
rection (j) in the Bloch sphere. Starting from these definitions, 
we construct the witness operator 



this entanglement witness coincides with the CHSH-based in- 
equality of Eq. ( |A3t , given that no assumption is made on the 
efficiency of the detection apparatus. 



W = C„K,x';<^') + C„(a;,x';^)+ 



(B7) 



In order to define the bounds on (W) satisfied by separable 
states, we consider a generic MMS-separable state described 
by the density matrix ff^^ ~ ^iPiPf ® pf ■ After detection 
losses on the multiphoton mode k^, such state evolves into 



1 ^ , which gives 



+ {AY{B%-{AY{B)];) 



(B8) 



where 



(i?>;=Tr 



o 



(qx,x;'7)Pj 



riB 



(i)'=Tr[a^(^)pf 
(i')''=Tr [a^ 



pf] 



(B9) 



As all these terms satisfy 
{i, we get 

IW^""! < 2, 



< 1 with X 



(BIO) 



which is the desired witness condition. We conclude by dis- 
cussing the features of this inequality. On one side, we note 
that the derivation of this bound is performed under the as- 
sumption that the state is measured with efficiency r/. Hence, 
such a witness operator permits us to demonstrate the entan- 
glement before detection losses. On the other side, no assump- 
tion is necessary on the MMS source due to the generality of 
the derived criterion. Finally, we note that for the case rj = 1 

I 



Appendix C: Correlator for the CHSH-based test in ideal 
conditions 



In this section we report the full calculation of the correlator 
C {X^ , , x; 0) reported in the main letter We begin with the 
two-mode correlation Q, defined as 



Q {a^ , ax_L , X; = (0) » ("x) 



(CI) 

This operator corresponds to the measurement of the general- 
ized parity operator on both polarization modes {tt^, tt^^ } of 
the macropart of our state. The average Q{ay^, , 0) = 
ab{'^^\Q\^^)ab is related to the correlator of the CHSH- 
based inequality by 



(C2) 



This expression holds by considering the closure relation 



{a^^ ) = l^j^ , which in turn comes from the 



normalization of the Wigner function. 

1. Two-mode correlator 

We now calculate the two-mode correlator 
Q{ax,ct^^,X;(l)). Let us recall the expression of the 
micro-macro state under investigation: 



|1>^)b) 



(C3) 



where the state has been expressed in a generic equatorial po- 
larization basis {7r0,7r0^}. The value of Q(a^, a^^ , x; 0) 
is obtained by exploiting the relation between the two-mode 
n^(a^) (E) Il^^{a^^) operator and the two-mode Wigner 

function B ($ I (a^ ) (8) (a^^ ) I s = ^ (a^ , ) . 
We get 



2(ax:"x_L>X; 



(C4) 



Here, WF and WF stand for the two-mode Wigner functions of amplified 



and 



evaluated at the rotated phase-space variables {a^, a^^ }. The correlator Q [a^, a^^ , x; ' 
expression of the Wigner functions j32!| (where S = sinh g and C — cosh g) 



single-photon states respectively, 
)) is then derived starting from the 



<J«x>«xJ = ^{4| 



|^ + |a^P)(l+2S=)+2CSRo(a2 



<K>«xJ--j{4| 



P(l + 2S'2) + 2CS'Re(a^^e*'^)] -l}e"^^ 
^(1 + 2S^) + 2CS Re{ay*)\ - l) e-2[(l".J^+l".l^)(i+2S^)+2CSRe(4^< 



(C5) 
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by rotating the polarization of the phase-space variables 

{a^,Q;0^} as 



= e*(x cos(x - 0) - i^x^ siii(x - 0)], 

= e^(x-0)/2[^^^ cos(x - 0) - siii(x - 0)]. 



(C6) 



Finally, we replace the complex phase-space variables with 
the real quadratures {X^^ P^^X^^ , P^^ ) and obtain the full 
expression for Q{X^^ P^^X^^ , P^^ , x; </>)• However, this is 
too lengthy and rather uninformative and will not be reported 
here. 



2. Single-mode correlator 

We now calculate the single mode correlator 
C(Xy^, P^,X;4>)- The choice of this measurement oper- 
ator allows us to capture the nonlocal features of the MMS 
state generated by amplification of an entangled photon pair 
To evaluate this quantity we exploit Eq. iC2i . 



C(Xx,Px,x;0)=: 



dnQ{X^,P^,X^^,P^^,x; 



(C7) 

where the integral in cPa^^ has been replaced by the integral 
in the quadrature variables dVl—dX^^dP^^. After straight- 
forward algebra, we obtain the following expression for the 
correlator 



cos[2(x - </>)]e-2(^ 



(C8) 



where {X^,P^} define a set of rotated variables X-^ = 

XxCos(x/2) - PxSin(x/2), = + 

PxCos(x/2). The maximum of such a correlation opera- 
tor is obtained at the origin of the phase-space and reads 

C(0,0,x;<^) = cos[2(x-0)]. 



Appendix D: Correlator for the CHSH-based test under 
detection losses and nonunitary injection efficiency 

Here we report in details the calculation of the correlator 
Cp ^, when detection losses and a nonunitary injection effi- 
ciency are taken into account. These two effects represent the 
two main issues for an experimental observation of entangle- 
ment in a MMS. 

The model for the effect of losses at the detection stage 
is performed by inserting a beam-splitter of transmittivity rj 
along the transmission path of the field on mode k^. The 
other port of this beam-splitter is injected with a vacuum 
state, thus introducing vacuum-noise fluctuations in the sys- 
tem. Here we demonstrate that the correlator C,^ in presence 
of detection losses 7] can be evaluated as the convolution of 



the lossless correlator C with a Gaussian function of the form: 

c„(x„p„x;</') =^^(1^ / / dKdPx c(x;,Pi,x;0; 

xe I J. 

(Dl) 

We begin by writing the density matrix p* of the micro- 
macro state after losses occur at the detection stage 

pr-^{l0)A(</'l®/:[l<i>'^-)B($^-|] 
+ |0±>a(0±|®/:[|$^>b($1] 

~\<j>j_)Amc[\<P^)B{<^^^\]}, 



(D2) 



where £[■] is the map that describes the action of detection 
losses. The evaluation of the correlation operator Q on this 
density matrix leads to 

2^("x' "x-L ,X;4') = Y ['^v,4>± ("x' "xi)-^^^0("x' "xi )] 

(D3) 

where W^^ and W^^^ are the Wigner functions of the 
macrostates 1$*^) and 1$*^^) after losses. The action of de- 
tection losses in the phase-space can be written in the form of 
a Gaussian convolution 0311 



{X,P) = j j dX'dP' W{X,P)Kr,{X,P,X',P'), 



where K,{X, P, X', P')^^ exp{-2[ <^-v^f '^^ 



(D4) 



— —]}■ The correlator is obtained from as 
2 

TT 



dnQ^{x^,p^,x^^,p^^,x;^) 

(D5) 

By writing explicitly the Wigner function after losses as a 
Gaussian convolution we obtain 

C,(x^,P^,x;</>) = ^ 1 1 dX^^dP^^i{x'^,P'^) 

(D6) 

where 



I(Xi,Pi) = j j dX'^JP'^^Q{X'^,P'^,X'^^,P^^,x-A) 
dX^^ dP^^ K,^{X^j_,P^^,X'-^^,P!^^). 

. (D7) 

By changing the integration variables as X^^^ — > X^^^ = 

plicit function 

I(X;,Pi) = 



P ^ P = 

^ X± ^ ^ XI- 



we have the ex- 



2\J\ 



7r(l - 77) 

dKJPLQiK'Px^K.^Px.^^'<f^)^ 

(D8) 
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where |J| = 1 — rj. Eq. ( IDll i is found by integrating 
over dX^^dP^^, using Eq. Q to have I{X'^,P'^) = 
C(X^, P^, x; ^) and replacing this in Eq. (ID6b . 

We now proceed with the explicit calculation of Eq. ( IDll i. 
As a first step, we rotate the quadratures and the 

integration variables {X'^, P^) as 



Px=A'x sin(x/2)+Px cos(x/2) 



(D9) 



with X = (X, X') and V = (P, P') and the convention that 
only primed (unprimed) variables are involved in the equa- 
tions above. The correlator C,^ can be then expressed as a 
function of the rotated variables. After replacing the expres- 
sion of K,j in the correlator C^, it is matter of some straight- 
forward (although tedious) algebra to find that 



cos[2(x- 0)16 



■ M 



'-FT 



1 - 



(l-7^)(l + 27?n) 
1 + 4r/(l — 



27, 



4?7(1 — T])n 



with M = 7?e^f + (1 - 77) and M = r\e "^^ 



(DIO) 
(1 - 77). This 



expression is maximized at the origin of the phase space, read- 
ing 

C^(X^,P^,X;<t>) = cos[2(x- 0)]/:(77,g) (DID 
where the loss function £(77, g) has the form: 



77 + 277(1 — 77)71 

(1 + 477(1 - 7])n)^/'^' 



In typical experimental conditions, the injection of the sin- 
gle photon of the entangled pair \iIj~)ab into the OPA occurs 
with an efficiency p < 1 because of the imperfect matching 
between the optical modes of the amplifier and the single- 
photon one. Such nonideality can be modeled by allowing 
for a probability p of correct single-photon injection and a 
complementary probability (1 — p) that just vacuum state is 
injected in the amplifier and no correlations between the two 
output modes are set. This modifies the density matrix of the 
output modes as 



pf =p\i^-)AB{^-\ + {l-p)^-(^\0)B{Ol (D13) 



where Ia = \H)a{H\ + \V)a{V\ is a completely mixed 
single-photon polarization state, and |0)b(0| is the vacuum 
state. The bipartite state after the ampHfication process then 
reads 



Pp 



:p|*-)ab(*-| + (1-p)^®(c/opa|0)b 



\U. 



t 

OPA 



(D14) 



We can now proceed with the calculation of C{a^, x; 0) as 
C (a;,, X; 0) = P AS |<T^ (0) ® (a^ , x) I > AB 



-(l-p)Tr 



Uopa\0)b{0\U, 



(D15) 

As the second term factorizes (due to the lack of quantum cor- 



relations) and Tr 



^^^(0) 



— 0, such contribution is null. 



Therefore, under nonideal injection efficiency, the correlator 
is related to the ideal one according to Cp(X^, P^, x; 0) = 
p C {X^, P^, x; 0)- This result can be extended to the case of 
nonunitary detection efficiency, leading to 



(DI6) 



Appendix E: Correlator for the entanglement witness after 
detection losses and nonunitary injection efficiency 

Here we sketch the steps needed for the calculation of the 
correlator Cp ,j entering the entanglement test based on the 
witness operator of Eq. jBlj under losses and nonideal photon 
injection. By using arguments similar to those put forward in 
the previous sections, we have 



C; (a^ , x; = i {Tr [£ [I ) B (1>*- 1] 6^,^, (.a^,X;v) 



- Tr 



(El) 

(D12) where £[■] is the map describing the lossy process. 

We focus on the case ?7 > ^. By exploiting re- 



sults that have been previously obtained here, we have 
We now exploit the chain of relations 



- J c?^ax-Ly (W7^(ax."xi) - ^^^(ax.axi)) 
^ J (fa-^^Q^{a^,a^^,x](t)) = C,,(a^, x; </>) 



so as to get (a^ , x; 0; = ^C,f ^ ("x ' X; ' 
ogous procedure, we obtain 



(E2) 

With an anal- 



2C,^(ax:X;'/') if??<5- ^^^^ 



We can further generalize this result so as to take into account 
the effect of a nonunitary injection efficiency and finally get 
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